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Abstract 

Relativistic prescriptions relating Sachs form factors to nucleon charge and 
magnetization densities are used to fit recent data for both the proton and the 
neutron. The analysis uses expansions in complete radial bases to minimize 
model dependence and to estimate the uncertainties in radial densities due to 
limitation of the range of momentum transfer. We find that the proton charge 
distribution, fitted to recent recoil-polarization data displaying an almost lin- 
ear decrease in Gep/Gmp for ^ I (GeV/c)^, is significantly broader than 
its magnetization density. We also find that the magnetization density is 
broader for the neutron than the proton. The neutron charge form factor is 
consistent with the Galster parametrization over the available range of Q^, 
but the relativistic inversion produces a softer radial density. Discrete ambi- 
guities in the inversion method are analyzed in detail. The method of Mitra 
and Kumari ensures compatibility with pQCD at large and is most useful 
for extrapolating form factors. Although a recent observation that QF2p/Fip 
is approximately constant for 2 < < 6 (GeV/c)^ appears to be inconsistent 
with the scaling expected from quark helicity conservation, our analysis 
fits these data while remaining consistent with pQCD for large Q^. 
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I. INTRODUCTION 



A fundamental test of the QCD confinement meclianism, either on the lattice or in 
models, is the electromagnetic structure of the nucleon. This electromagnetic structure is 
reflected by the electric and magnetic form factors, GEiQ"^) and Gm{Q'^), measured by elec- 
tron scattering with spacelike invariant momentum transfer Q. At low Q one can interpret 
these form factors as Fourier transforms of the nucleon charge and magnetization densities 
but the relationship between form factor and density is complicated by recoil as Q 
increases. Although models of nucleon structure can often calculate the form factor directly, 
it is still desirable to relate form factors to spatial densities because our intuition tends to 
be grounded more firmly in space than momentum transfer. In this paper we fit charge 
and magnetization densities to recent nucleon form factor data using a prescription that ac- 
counts for nucleon recoil and Lorentz contraction and is compatible with perturbative QCD 
(pQCD) scaling at large Q"^. 

Early experiments with modest suggested that 

Gmp Gmu ^ /.x 
Gep ~ ~ ~ Gd (1) 

where GoiQ^) = (1 + Q^/A^)-^ j\2 ^ q ^i (GeV/c)^ is known as the dipole form 
factor P,^. However, the naive Fourier transform of the dipole form factor produces an 
exponential density with an unphysical cusp at the origin. Similarly, data for Geu at low 
can be described by the Galster parametrization ^ 

At 

GEnm ^ -f^uGoiQ')-—^ (2) 

i + IDT 

where A and B are constants and r = {Q/2m)'^, but direct Fourier transform of this form 
factor also produces a cusp at the origin. Using a relativistic inversion formula that ac- 
counts for the Lorentz contraction along the momentum transfer, Licht and Pagnamenta 
0] obtained a reasonable fit to proton form factors using a Gaussian density with a more 
realistic behavior in the interior. Ji obtained similar results using a relativistic soliton 
model and we used this model in Ref. to fit data for Sachs form factors. These models 
offer plausible radial densities, but are not compatible with pQCD scaling unless one im- 
poses somewhat awkward restrictions upon the Fourier transform, as discussed in Sec. [11 C 



Fortunately, a variant proposed by Mitra and Kumari complies with pQCD scaling au- 
tomatically, without need of such constraints. We use this method, described as relativistic 
inversion, to extract nucleon charge and magnetization densities from data for Sachs form 
factors. Our fitting procedure minimizes model dependence by employing linear expansions 
in complete radial bases, such as Fourier-Bessel or Laguerre-Gaussian expansions, and es- 
timates uncertainties arising from the limitation of experimental data to a finite range of 
momentum transfer using methods originally developed to analyze electron scattering by 
nuclei. Such an analysis produces a good fit to form factor data using a radial density 
whose error band reflects both the statistical quality of the data and its limited coverage of 
momentum transfer. Differences between densities obtained using several variations of the 
inversion formula are described as discrete ambiguities and are analyzed in detail herein. 
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Data for Gmp and Gmu with > 1 (GeV/c)^ show significant departures from the 
simple dipole parametrization, but the extraction of Gep using the traditional Rosenbluth 
method [jlO[ becomes increasingly difficult as increases because the dominance of the Gmp 



contribution to the cross section increases with . Consequently, there are large statistical 
uncertainties in Rosenbluth data for Gep at > 1 (GeV/c)^ and the discrepancies between 
comparable experiments suggests that systematic errors in the Rosenbluth analysis are often 
underestimated |TII]. More recently, recoil polarization has been used to measure the ratio 



Qp = Gep/Gmp directly, without need of Rosenbluth separation. In this technique, the 
components of the nucleon polarization P' after scattering by a polarized electron beam 
are measured along the momentum transfer direction, denoted by 5, and in the x direction 



transverse to z in the scattering plane. The form factor ratio is then obtained using |12.13 



PL 2e 



X 



r 1 + e 



9 (3) 



where e = (1 + (1 + r)2 tan^ 9e/2)~^ is the transverse polarization of the virtual photon for an 
electron scattering angle Of.- For the proton, both components can be measured simultane- 
ously using a polarimeter in the focal plane of a magnetic spectrometer, thereby minimizing 
systematic uncertainties due to beam polarization, analyzing power, and kinematic param- 
eters. The systematic uncertainty due to precession of the proton spin in the magnetic 
spectrometer is usually much smaller than the systematic uncertainties in combining the 
absolute cross sections obtained with different kinematical conditions and acceptances that 
are needed for the Rosenbluth method. Recent data using the recoil polarization technique 
T^ , p!5|JT6| have shown a dramatic, almost linear, decrease in Gep/Gmp for > 1 (GeV/c)^. 



It was suggested that those results demonstrate that the proton charge is distributed over a 
larger volume than its magnetization, but radial densities were not obtained. Our analysis 
confirms that interpretation quantitatively. 

Similar techniques can be used to obtain the neutron form factor ratio, Qn = GEnlGum 
using either target or recoil polarization, but in the absence of a target with free neutrons 
one must employ quasifree scattering from a neutron bound in a light nucleus. Detection of 
a recoil neutron with quasifree kinematics and small missing momentum tends to minimize 
uncertainties due to nuclear structure and final state interactions ITT]. Although considerable 



care is still needed at low Q^, polarization methods offer smaller systematic errors and less 
model dependence than traditional Rosenbluth analyses of elastic scattering or quasifree 
knockout. We extracted the neutron charge density from recent polarization data for for 
< 1.6 (GeV/c)^ using the relativistic inversion method. Although the form factor data 
remain consistent with the Galster parametrization over this range of momentum transfer, 
the charge density obtained by relativistic inversion is considerably softer than that from 
nonrelativistic inversion of the Galster form factor and does not feature a cusp at the origin. 
Over the next several years, extending of the experimental range of momentum transfer 
should substantially reduce the uncertainty in the interior density. 

The model is presented in Sec. 0, the analysis procedures in Sec. |T|, and principal results 
in Sec. fV[ In Sec. |V| we compare our results to another analyses, discuss the extrapolation 
to higher and the role of discrete ambiguities in fitted densities. We also form combi- 
nations of neutron and proton charge densities that in the naive quark model represent the 



3 



distribution of up and down quarks in the proton. Finally, our conclusions are summarized 
in Sec. IVI, 



II. MODEL 

A. Sachs Form Factors 

Matrix elements of the nucleon electromagnetic current operator take the form 

{N{j)\ s')\J^'\N{j), s)) = u{p', s')eT^'u{p, s) (4) 

where u is a Dirac spinor, p,p' are initial and final momenta, q = p — p' is the momentum 
transfer, s, s' are spin four-vectors, and where the vertex function 

= F,{Q')Y + ^F,{Q')'^ (5) 

2m 

features Dirac and Pauli form factors Fi and F2 that depend upon the nucleon structure. 
Here e is the elementary charge, m is the nucleon mass, k, is the anomalous part of the mag- 
netic moment, and 7^ and a'^'^ are the usual Dirac matrices {e.g., |18|). The interpretation 
of these form factors appears simplest in the nucleon Breit frame where the energy trans- 
fer vanishes. In this frame the nucleon approaches with initial momentum —qB/2, receives 
three-momentum transfer qb, and leaves with final momentum 5*^/2. Thus, the nucleon 
Breit frame momentum is defined by Qb ~ Q"^ ~ + where {uj,q) is the momentum 

transfer in the laboratory, = — uj"^ is the spacelike invariant four-momentum transfer, 
and T = Q'^ /Am?. In the Breit frame for a particular value of Q^, the current separates into 
electric and magnetic contributions 

u{p\ s')T^u{p, s) = xi {ge + '-^^Gm^ Xs (6) 

where Xs is a two-component Pauli spinor and where the Sachs form factors are given by 

Ge = Fi- TKF2 (7a) 
Gm = Fi + ^F^ (7b) 

The similarity of Eq. (^ to the classical current density 

^^^ = (epr,/^-xvpr) (8) 

suggests an identification of charge and magnetization densities 

2 roo 



') POO 

pT{r) = - dQQ'UQr)GE{Q') (9a) 

TT JO 

/^Pr(0 = - / dQ Q'UQr)GMm (9b) 

TT Jo 



where /i = 1 -|- k is the appropriate static magnetic moment (either /ip or //„) relative to the 
nuclear magneton. However, this naive inversion procedure is described as nonrelativistic 
(NR) because it ignores the variation of the Breit frame with Q^. 
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B. Intrinsic Form Factors 



Let Pch{r) and Pm{f) represent spherical charge and magnetization densities in the nu- 
cleon rest frame. It is convenient to normahze these intrinsic densities according to 



/ dr r^pchir) = Z (10a) 
Jo 



oo 







dr r>„(r) = 1 (10b) 



where Z = 0, 1 is the nucleon charge. Fourier-Bessel transforms of the intrinsic densities are 
defined by 



oo 



p{k) = / dr r jolkr) p{r) (11) 







where k is the spatial frequency (or wave number). We describe p{k) as an intrinsic form 
factor. If one knew how to obtain p{k) from data for the appropriate Sachs form factor, the 
intrinsic density could be obtained simply by inverting the Fourier transform, such that 

2 

p{r) = - dk k'Ukr)p{k) (12) 



TT JO 

The naive nonrelativistic inversion method assumes that k Q and p{Q) — > G{Q'^) 
where G{Q^) is the appropriate Sachs form factor. However, this inversion procedure pro- 
duces unsatisfactory results for the common dipole and Galster parametrizations — the 
corresponding radial densities have unphysical cusps at the origin and rather hard cores. 
For example, the naive Fourier transform of the dipole form factor produces an exponential 
density. (Although it appears much more complicated, the Galster density can also be ob- 
tained in closed form and displays similarly unrealistic behavior near the origin.) Licht and 
Pagnamenta [0 attributed these failures of nonrelativistic inversion to the replacement of 



the intrinsic spatial frequency k with the momentum transfer Q and demonstrated that by 
applying a boost from the Breit frame with momentum qs = Q to the rest frame, inversion 
of the dipole form factor using a reduced spatial frequency 

e ^ ^ (13) 

softens the density. In fact, a good fit to the data for Gep could then be obtained using a 
Gaussian density typical of quark models. 

Unfortunately, unique relativistic relationships between the Sachs form factors measured 
by electron scattering at finite and the static charge and magnetization densities in the 
nucleon rest frame do not exist. The basic problem is that electron scattering measures tran- 
sition matrix elements between states of a composite system that have different momenta 
and the transition densities between such states are different from the static densities in 
the rest frame. Furthermore, the boost operator for a composite system depends upon the 
interactions among its constituents. Nevertheless, a wide variety of models have employed 
similar relativistic prescriptions to relate elastic form factors to ground-state densities. The 
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first proposal was made by Licht and Pagnamenta |T9[ using a cluster model and a kine- 
matic boost that neglects interactions. The transition form factors were evaluated using 
the impulse approximation and neglecting relative motion. Mitra and Kumari obtained 
similar results using a kinematic transformation that is more symmetric between initial and 
final states and can be applied to inelastic scattering also. Ji also obtained similar results 
using a relativistic Skyrmion model based upon a Lorentz invariant Lagrangian density for 
which the classical soliton solution can be evaluated in any frame. Quantum fluctuations 
were then evaluated after the boost. Although an approximation is still required to evalu- 
ate the transition form factors, it was argued that this approximation is best in the Breit 



frame. Holzwarth |2y] extended the soliton model to the timelike regime and analyzed the 
super convergence relations needed to obtain spectral functions. 
Each of these prescriptions can be represented in the form 

~Pch{k) = GE{Q^){l + rf'' (14a) 
^^pm{k) = GM{Q^){l + rf" (14b) 

where G{Q'^) is one of the four Sachs form factors, k is the intrinsic spatial frequency given 
by Eq. (|T3p, and A is a model-dependent constant. The most important relativistic effect is 
Lorentz contraction of spatial distributions in the Breit frame and the corresponding increase 
of spatial frequency represented by the factor of (1 + r) in Eq. (|l^). A measurement with 
Breit-frame momentum transfer Qb = Q probes a reduced spatial frequency k in the rest 
frame. The Sachs form factor for a large invariant momentum transfer is determined by 
a much smaller spatial frequency k"^ = Q^/(l -l- t) and thus declines much less rapidly with 
respect to than the Fourier transform of the density declines with respect to k"^. In fact, 
the accessible spatial frequency is limited to k < 2m such that the asymptotic Sachs form 
factors in the limit Q"^ ^ oo are determined by the Fourier transform of intrinsic densities in 
the immediate vicinity of the limiting frequency km = 2m. In this model, no information can 
be obtained beyond the limiting frequency determined by the nucleon Compton wavelength. 
This limitation can be understood as a consequence of relativistic position fluctuations, 
known as of zitterbewegung, that smooth out radial variations on scales smaller than the 
Compton wavelength. 

Ji [0 derived A^; = for electric and Am = 1 for magnetic form factors in the soliton 
model and attributed the difference between A^; and Am to the Lorentz transformation 
properties of scalar and vector densities. The same choices were employed by Holzwarth 
pO| , |2T[| . On the other hand, Licht and Pagnamenta [|19[ obtained A^; = Am = 1 using the 



cluster model, but Mitra and Kumari ^ found that a more symmetric treatment of the 
kinematics gives values A^; = Am = 2 that automatically satisfy the perturbative QCD 
scaling relations at very large Q^. For most of the present analysis we will employ the 
symmetric choice A^; = Am = 2 because fewer restrictions upon the behavior of p{k) are 
needed near km to ensure compatibility with the asymptotic behavior of Sachs form factors 
expected from dimensional scaling. We describe the variation of a fitted density with the 
choice of A as a discrete ambiguity. The effect of discrete ambiguities upon fitted densities 
will be examined in Sec. |V D| . 
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C. Asymptotic Behavior 



The present form factor model with A > suggests that the asymptotic behavior for 
Q ^ 2m is given by 



2Q2 

+ (4A(A + l)p{kj + (3 + 4X)k^p'{k^) + klp"{k^))^ + ...^ (15) 

where p' and p" are derivatives of the momentum-space density evaluated at the limiting 
frequency km = 2m. Evidently, noninteger values of A are incompatible with the perturbative 
QCD prediction [^,^ that G ~ Q~'^ aside from logarithmic corrections. Similarly, A > 2 
is excluded also, leaving just three choices. If we choose A = 2, then we need only require 
p{km) 7^ to obtain consistency with pQCD. Thus, the proposal by Mitra and Kumari 
of A^; = Am = 2 offers the most natural approach to the pQCD limit. If we choose A = 1, 
as recommended by Refs. |]7|j20[| for Gm or by Ref. |]19[ for both Ge and Gm, then we must 



require p{k) to have a node at km such that p{km) = and p'{km) 7^ 0. Finally, if we choose 
A = as recommended by Refs. ||7|,pO| for Ge, then we must impose the somewhat unnatural 
constraints p{km) = p'{km) = with p"{km) 7^ 0. Thus, it appears that the usefulness of the 
chiral soliton model is limited to <^ 4m^ and in order to fit data for larger with that 
model Holzwarth found it necessary to artificially increase the soliton mass 

Although the intrinsic form factors p{k) obtained using either dipole or Galster functions 
for G{Q^) in Eq. (0) are compatible with the pQCD constraints upon p{km), neither can 
be inverted using Eq. (JT^) with A = because p{k) ~ A;^)(— /c^/^)^^ for k ^ 00. The 
inversion integrals for these functions converge well for A = 2, slowly for A = 1, and diverge 
for A = 0. Recognizing that pQCD favors A = 2, we expect p{k) to have an asymptotic 
k~'^ behavior with an amplitude determined by the nucleon-antinucleon annihilation process 
NN — i> e~e~^ at threshold. Similarly, the behavior of p{k) for /c — 00 should be determined 
by the electromagnetic annihilation data for ^ — 4m^. However, we have not attempted 
to incorporate electromagnetic data for timelike in the present analysis because it is not 
clear that the prescription for intrinsic form factors should apply to that regime. A more 
general analysis of the analytic structure of the form factors can be made using dispersion 
theory pi| , pB| , pB| , but that approach does not consider the densities that are the subject of 
the present analysis. Nor do we consider here the modifications of pQCD scaling due to 
logarithmic running of the strong coupling |p7 |. 



With the exception of Gmp, the available data constrain p{k) very little near km because 
the ratio 

A ~ 1 _ 3lL (16) 

km 2Q 

approaches unity relatively slowly as increases. Thus, we will find that the choice of 
A has very little effect upon the fit to data for Sachs form factors, but does have a strong 
effect upon the extrapolation of fitted form factors beyond the measured range of Q^. By 
incorporating pQCD scaling in its basic parametrization, the choice A = 2 limits the range 
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of variation available to extrapolated form factors. Conversely, without explicit enforcement 
of pQCD by means of somewhat artificial constraints upon p{k) near km, fits with A < 2 
permit much wider latitude at large Q^. The data for Gmp for > 20 (GeV/c)^ exhibit 
scaling and automatically enforce the appropriate constraints upon p{kjn), but data available 
for the other three electromagnetic form factors do not. Consequently, one could impose 
constraints upon p{km) with little effect upon the fits in the measured range of Q^. However, 
we chose not to employ constraints of this kind, which seem rather artificial, and to permit 
fits with A < 2 the greatest possible latitude. 

We use A^; = Am = 2 for most of the present work, but will discuss the consequences 
of the discrete ambiguity in Sec. |V h)| . Note that our previous work 0, motivated by the 
soliton model, used A^; = and Am = 1- 



D. Moments 

It is customary to describe the low behavior of a form factor in terms of a transition 
radius obtained from integral moments of the underlying density, but care must be taken 
with the relativistic relationship between a Sachs form factor and its intrinsic density. We 
define integral moments by 

roo 

M„ = / dr r2+"p(r) (17) 
Jo 

where a is an even integer. For a charge density these moments are related to the electric 
form factor by 

Mo = G(0) (18a) 



Ma = -6 ' 



G{0) (18b) 



2m2 



while for magnetization we divide by the magnetic moment. Thus, one expects Mq = 
Z for charge densities and Mq = 1 for magnetization densities. Notice that the lowest 
nonvanishing moment is free of discrete ambiguities, but that higher moments depend upon 
A. For example, the mean-square neutron charge radius reduces to 



. 2^ . dG{Q') 
{r )n = -6 



(19) 



dQ^ 

because the charge vanishes, while the proton radius retains a small dependence upon A 

= {r\„ - 1^ (20) 

due to the discrete ambiguity in the intrinsic density. This term, equal to 0.066A in units 
of fm^, appears to be similar to the famous Foldy contribution to the neutron charge radius 
p8| but has a different origin because it does not depend upon the anomalous magnetic 
moment. 
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When Mo 7^ it is useful to distinguish between a radius parameter 



based upon the initial logarithmic derivative of a Sachs form factor and the rms radius 

^rr.--(e-^f (22) 

of the corresponding density obtained for specified A. Thus, ^ is a model-independent 
property of the form factor data while (r^)^/^ is subject to a discrete ambiguity. These 
radii agree for A = 0, but (r^)^/^ is smaller than ^ for larger A due to the zitterbewegung 
correction. 

Accurate calculations for many phenomena in atomic physics, such as the Lamb shift, 
require corrections for the finite size of nucleons. Although it might appear that the nucleon 
size should be determined by an integral moment of a nucleon density through Eq. (|1^ , the 
static radial density is not directly measurable by electron scattering. The discrete ambigu- 
ity between the initial slope of the Sachs form factor and its associated transition density 
reflects the model-dependence of the relativistic inversion procedure arising from the treat- 
ment of zitterbewegung. By convention, QED theorists have decided to identify the radius 
with the initial slope of the Sachs form factor and to treat recoil, vacuum polarization, zit- 
terbewegung, and other effects as separate corrections {e.g. |2^j3^). To distinguish between 
various determinations of nucleon size, we describe the model-independent quantity ^ as the 
Sachs radius and the model-dependent rms radius obtained from moments of a fitted radial 
density as an intrinsic radius. When necessary, these radii are further qualified as either 
charge or magnetic and for intrinsic radii by the value of A. The Sachs charge radius is 
usually the most appropriate for QED applications. 



E. Example: Gaussian density 

It is instructive to consider the Sachs form factor that would be produced by a simple 
Gaussian density 

p(r) = exp (-(r/6)2) ~p{k) = exp {-{kb/2f) (23) 

that is typical of quark models. The form factor obtained using Eq. (14) with A = 2 is 



compared with the familiar dipole form factor in Fig. |1] for several choices of b; note with 
b = 0.556 fm the Gaussian parametrization has the same rms radius as the dipole form factor. 
These curves display the same general features as the data for Gep, Gmp and Gmu- for low 
the form factor is close to the dipole form while for large one finds an asymptotic 
limit for G/Gc that depends sensitively upon b but is less than unity for reasonable values. 
The greatest sensitivity to the shape of the density is found in a transition region for that 
ranges from several tenths to several (GeV/c)^, depending upon b. Thus, data with similar 
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general features can be fit by modulating a basic Gaussian with an even polynomial, where 
the polynomial degree can be minimized by an optimal choice of b. For Geu one need only 
require the polynomial part of p{k) to begin with fc^ to ensure that the net charge vanishes. 
Expansions of this form are no more complicated than other parametrizations in common 
use, but are free of unphysical cusps at origin. 



III. ANALYSIS PROCEDURES 

A. Linear expansions 

To extract radial densities from the nucleon form factor data we employ techniques 
originally developed for fitting radial distributions to data for scattering of electrons or 
protons from nuclei pT| , |32| , |33| . Simple models with a small number of parameters do not 
offer sufficient fiexibility to provide a realistic estimate of the uncertainty in a radial density. 
Rather, we employ linear expansions in complete sets of basis functions that are capable 
of describing any plausible radial distribution without strong a priori constraints upon its 
shape. Such expansions permit one to estimate the uncertainties in the fitted density due 
to both the statistical quality of the data and the inevitable limitation of experimental data 
to a frequency range, k < kmax, where 

k - ^"^"^ (OA) 
+ ^^^^ 

is the maximum spatial frequency sampled by experimental data limited to Q < Qmax- The 
uncertainty due to limitation of k is known as incompleteness error. 
A radial density can be represented as an expansion of the form 

PW=E«n/n(0 (25) 

n 

where the basis functions /n(r) are drawn from any convenient complete set. The corre- 
sponding Fourier transform then takes the form 

m=T.^nfnik) (26) 

n 

where 

Uk)= / drr^3,{kr)Ur) (27) 

<J 

represent basis functions in momentum space. The expansion coefficients, fitted to 

form factor data subject to several minimally restrictive constraints to be discussed shortly. 
Analyses of this type are often described as model independent because a complete basis 
can reproduce any physically reasonable density; if a sufficient number of terms are included 
in the fitting procedure the dependence of the fitted density upon the assumptions of the 
model is minimized. By contrast, simple parametrizations like the Galster model severely 
constrain the shape of the fitted density. 
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We consider two bases that have been found useful in the analysis of electron or proton 
scattering data. The present discussion is limited to monopole densities, but generalizations 
to higher angular momenta are discussed in Refs. p3| , p^ . 

The Fourier-Bessel expansion (FBE) employs basis functions of the form 

/n(r) = Uhnr)Q{Rma. ' r) (28a) 

fM = ^^^^J0{kRma.) (28b) 

where 6 is the unit step function, Rmax is the expansion radius, and kn = mr/Rmax- One 
advantage of the FBE is that the contribution of each term to the form factor is concentrated 
around its characteristic frequency /c„ so that a coefficient a„ is largely determined by data 
with k ~ kn- The larger the expansion radius Rmax, the smaller the spacing between 
successive fc„ and the greater the sensitivity one has to variations in the form factor. One 
should choose Rmax to be several times the root-mean-square radius but not so large that an 
excessive number of terms is needed to span the experimental range of momentum transfer. 
Terms with kn > kmax provide an estimate of the incompleteness error. We chose Rmax = 4.0 
fm, but the results are insensitive to its exact value. However, a disadvantage of the FBE 
is that a relatively large number of terms is often needed to accurately represent a typical 
confined density. 

Alternatively, the Laguerre-Gaussian expansion (LGE) employs basis functions of the 
form 

fn{r) = e-^'L]!\2x^) (29a) 
= ^h\-re-y'L]/\2y') (29b) 

where x = r/b, y = kb/2, and is a generalized Laguerre polynomial. A significant advan- 
tage of the LGE is that the number of terms needed to provide a reasonable approximation 
to the density can be minimized by choosing b in accordance with the natural radial scale. 
We chose b = 0.556 fm such that the mean-square radius of the Gaussian factor is consistent 
with that of the common dipole parametrization of Sachs form factors. We then find that 
the magnitude of decreases rapidly with n, but the quality of the fit and the shape of 
the density are actually independent of b over a wide range. However, a disadvantage of the 
LGE is that the basis functions are not localized in momentum space so that the coefficients 
tend to be correlated more strongly than for the FBE. 

B. Constraints 

The expansion coefficients are obtained by minimizing 

where is the fitted value of a quantity yi with uncertainty 6yi. In addition to experimental 
data, the set yi generally includes pseudodata used to enforce constraints and to estimate 
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the incompleteness error associated with the hmitation of experimental data to a finite range 
of momentum transfer. 

The absence of data for very large requires some constraint upon the behavior of 
p{k) for kmax < k < km- Furthermore, inversion of the Fourier transform also requires an 
assumption about the experimentally inaccessible region k > km- On quite general grounds 
one expects the asymptotic form factor for a confined system to decrease more rapidly than 
k~^ in particular, this condition ensures that there will be no cusp at the orgin. In fact. 



our results show that the intrinsic form factor p(k) is well approximated by a Gaussian for 
large k. Therefore, we will assume that p{k) for k > k^, is bounded by a envelope and 
use the fiexibility afforded by that envelope to estimate the incompleteness error due to the 
limitation of experimental information to the range k < kmax- Although some restriction 
is needed to stabilize the fits, the fc""^ envelope probably overestimates the uncertainties in 
unmeasured form factors and their effect upon uncertainties in fitted densities; nevertheless, 
we prefer to employ minimally restrictive constraints so that those densities will have the 
best possible model independence. 



More detailed discussions of the method may be found in Refs. pl|j3^ , |33[| , but the basic 
idea is to supplement the experimental data by pseudodata of the form p{ki) = ± 5p{ki) 
whose uncertainties are based upon a reasonable model of the asymptotic behavior of the 
form factor for ki > k^ax where kmax is the spatial frequency corresponding to the maximum 
measured Q^. Therefore, uncertainties in the form factor for k > kmax are based upon an 
envelope of the form 

Sp{k) = Piim{k) (31a) 

Plim{k) = \pikmax)\ l^^j (31b) 



where the factor of yl/S represents the variance of a uniform distribution of unit width. 
When using FBE the pseudodata are chosen at the characteristic frequencies fc„ = nn/Rmax 
with n > kmaxRmax/T^i while a uniform spacing of Afcj = 1.0 fm~^ was employed for LGE. 
The error band for a fitted density is computed from the covariance matrix for the fit and 
includes the incompleteness error. A detailed discussion of the decomposition of the density 
uncertainty into statistical and incompleteness errors may be found in Ref. . 



Recognizing that pQCD imposes an asymptotic limit of the form G oc upon the 
Sachs form factors, one might be tempted to employ pseudodata for G at large Q^- If one 
knew how to estimate the proportionality constant, this procedure could be used to regulate 
p{k) for kmax < k < km but would not be sufficient for construction of the radial density 
because inversion of the Fourier transform also requires information for the k > km region 
that is inaccessible to electron scattering. Although we expect p{k) to be small for k > 2m, 
we cannot simply set it to zero because an abrupt cut-off would introduce unreasonable 
density oscillations at very large radii. The present procedure estimates the uncertainty in 
the radial density arising from both the unmeasured and the unmeasurable ranges of spatial 
frequency. In this model, the minimum uncertainty in density is governed by the nucleon 
Compton wavelength and can be interpreted as an irreducible smearing by zitterbewegung. 
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Small but undesirable oscillations in fitted densities at large radii were suppressed using 



a tail bias based upon the method discussed in Ref . |3J] . We employed a tail function of the 
form t(r) oc e~^^, based upon the successful dipole parametrization for low Q^, and included 
in the fit a penalty function of the form 

,;.gf^!MziMV (32) 



i=l 



to suppress strong deviations from the tail function. The radial pseudodata were constucted 
on the grid rj = + ^Ar for i = l^N^ in the range r > r^- We chose A = 4.27 fm~^, 

= 2.0 fm, Ar = 0.2 fm, = 10, and w = 2, but the results are rather insensitive to 
these details. The tail bias improves the convergence of moments of the density but has 
little effect upon a fitted density in the region where it is large. 

The fitting procedure also permits constraints to be placed on integral moments of the 
radial density. We define fitted moments by 



oo 

2+a 



M^= dr r^+°p(r) (33) 







where here p(r) is the fitted density and include a penalty function of the form 

where is the measured value and 5Ma is its uncertainty. The constraint on the neutron 
charge was enforced by means of a pseudodatum Mq = ± 10^®. In addition, the atomic 
physics datum for M2 from Ref. was included in fits made to the neutron charge density. 
It is also useful to define a fitted transition radius R as 



R = VM2/M0 (35) 

for Gep, Gmp, or GMn- Thus, the fitted transition radius is correlated with the experimental 
normalization at low Q^. If the fitted Mq were constrained, the uncertainty in R would be 
artificially reduced. Therefore, no constraints were placed on Mq for Gsp, Gmp, or Gmu — 
given that those intrinsic densities were defined with unit normalization, the fitted values of 
Mq test the normalization of the experimental data. 



C. Data selection 

We tried to select the best available data in each range of Q^, with an emphasis upon 
recent data using recoil or target polarization wherever available. These selections are sum- 
marized in Table |. Although a thorough review of the data for nucleon electromagnetic form 
factors is beyond the scope of the present work, in this section we provide brief explanations 
for some of our selections and omissions. 

Gmp data were taken from the compilation of Hohler et al. [0] for Q"^ < 0.15 (GeV/c)^ 



and for larger Q from the compilation and re-analysis made by Brash et al. |36| to improve 
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the correction for the small contribution of Gep to cross section at large Q^. Values for Gep 
were obtained by multiplying the recent recoil polarization measurements of Gep/Gmp from 
MUM 



Refs. 



J8[ by the Brash parametrization of Gmp- Supplementary recoil polarization 
data from Ref. were omitted — those data are consistent with those selected but have 



larger statistical uncertainties. In addition, the cross section data for Gep from Refs. p9|j40 
were used at low Q^. The cross section data for Gep at > 1 (GeV/c)^ from Refs. ||4l|j42 
were omitted because, as shown in Ref. |T^, they are in significant disagreement both 



with the recoil polarization data and with each other presumably because the Rosenbluth 
technique becomes increasingly difficult as increases and the relative contribution of the 
electric form factor to the unpolarized cross section becomes quite small. Nevertheless, we 
eagerly await new results from a proposed improvement of the Rosenbluth method . 

The neutron magnetic contribution is large enough at high to employ quasielas- 
tic electron-deuteron scattering with subtraction of the proton contribution. Data of this 
type were obtained from Refs. [0,^. At low the model dependence of the quasielas- 
tic method becomes relatively large. Markowitz et al. |^6|] measured the quasifree neutron 
knockout cross section for the (i(e, e'n) reaction and calibrated the efficiency of the neutron 
detector using associated particle production in the deuteron photodisintegration reaction, 
(i(7,p?T,). The dependence upon the deuteron wave function can be reduced by analyzing 
the ratio between quasifree cross sections for neutron or proton knockout, indicated as the 
(i(e, e'N) reaction in Table |, with relatively small corrections made for meson-exchange cur- 
rents and final-state interactions (FSI). Bruins et al. calibrated their neutron detector 
using the ^(7, vr+n) reaction while Refs. |^8| , ^ , p0| employed elastic neutron-proton scat- 
tering. The associated-particle techniques permit calibration in situ but must correct the 
bremsstrahlung measurements for the contributions of three-body reactions from electropro- 
duction that lie outside the acceptance pT] , ^ |, whereas the p{n,p)n reaction is kinematically 
complete but requires calibration at a different facility and under different conditions than 
used for the reaction of interest. Unfortunately, these methods remain in substantial dis- 
agreement; because we are not convinced there is a compelling preference, we include the 
data from both methods in the present analysis. Finally, we also include recent low 
data from inclusive electron scattering from transversely polarized ^He that uses the 
Fadeev calculations of Refs. ||5^j55|l to correct for nuclear structure effects. These data are 



more consistent with the coincidence ratio method calibrated by elastic scattering than by 
associated-particle production. 

We do not use any elastic or quasielastic cross section data for Geu because the uncer- 
tainties arising from nuclear structure are prohibitively large. Polarization techniques offer 
a signal that is linear in Geu and with less model dependence. Nevetheless, at low it 
remains important to correct recoil polarization data for final-state interactions and target 
polarization for nuclear structure. Most of the data for deuterium targets have been ana- 
lyzed using the calculations of Arenhovel et al. 0,0] to correct for nuclear structure; in 
Table H we cite both the experimental paper and the subsequent analysis. The result from 
the first experiment of this type |Q has not been corrected, but the statistical uncertainty 
was large. The = 0.4 (GeV/c)^ data for ^He(e, e'n) have been analyzed using Fadeev cal- 
culations but Rohe et al. |^ argue that the plane-wave impulse approximation (PWIA) 
is adequate at larger Q^. We also include values for Geu extracted by Schiavilla and Sick 
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|60| from an analysis of the deutron quadrupole form factor obtained from tensor polar- 
ization measurements of elastic scattering |^,^,^,64|. Although the model uncertainties 
remain somewhat large, this analysis covers a larger range of and appears to be more 



accurate than the older analysis of deuteron elastic scattering by Platchkov et al. , which 
is omitted. 

Finally, the neutron mean-square charge radius is related to the neutron-electron scat- 
tering length, 6„e, by 



Q;m„c 



(36) 



Unfortunately, the measurements are rather difficult and most techniques require substantial 
corrections for effects for which there is often insufficient information. A recent review 
of these measurements has been made by Alexandrov [Q, who finds that most modern 
measurements cluster around two values. From measurements of the energy dependence for 
the transmission of thermal neutrons through liquid ^°^Pb, Kopecky et al. Ref. obtained 
bne = (-1.33 ±0.027 ±0.03) x 10"^ fm, corresponding to {r^)n = -0.115 ±0.003 fm^. This 
result agrees well with similar measurements by Koester et al. [37| for lead isotopes and ^''^Bi 



obtained bne = (-1.60 ± 
using neutron diffraction 



and with the results of Krohn and Ringo |^ using the angular distribution for neutron 
scattering by noble gases. Alternatively, Alexandrov et al. 
0.05) X 10~^ fm, corresponding to (r^)„ = —0.138 ± 0.004 fm 
from single crystals of ^^^W. This result is consistent with a bismuth transmission experiment 
that was also performed at Dubna, but disagrees by about 5 standard deviations from the 
Garching, Argonne, and Oak Ridge experiments. This discrepancy has been attributed to 
resonance corrections |^0[^ but remains controversial. The extracted scattering length is 
strongly correlated with the resonance correction. Leeb and Teichtmeister argue that 



the correction employed by Alexandrov et al. requires implausibly large contributions 
from negative energy levels. Alexandrov |7T|] argues that the energy-independent resonance 
correction should be fitted to the data and that in the absence of definitive information 
a negative contribution cannot be excluded. Futhermore, neutron diffraction from single 
crystals of ^'^^W provides a larger signal than the energy dependence of the total cross 
section, and is hence is less sensitive to this correction. Although this controversy has not 
yet been resolved satisfactorily, we decided to employ the most recent result from Oak Ridge, 
namely = —0.115 ± 0.003 fm^, as a datum in our fit of the neutron charge density and 
to omit the Dubna result. The sensitivity to this choice is discussed in Sec. [IV (J4| . 



IV. RESULTS 
A. Form Factors 

Fits to the form factor data are shown in Fig. |^ as bands that represent the uncertainties 
in the fitted form factors. These bands were computed using the covariance matrix. The fits 
shown in Fig. |^ employ the LGE parametrization with Xe = Xm = but the results using 
the FBE parametrization are practically indistinguishable. Nor do these fits depend upon 
the choices for b, Rmax, or details of the constraints. Fits using A < 2 are almost identical 
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within the ranges spanned by experimental data, but their error bands grow more rapidly 
at larger Q"^ . The rapidly decreasing dipole form factor is divided out to emphasize the 
deviations at large from this characteristic behavior. For Geu we also display a simple 
two-parameter fit using the Galster parametrization. 

The intrinsic form factors, obtained via the relativistic transformation prescribed by Eq. 
(0), are shown in Fig. ^ using A^; = Aa/ = 2. From these figures we observe that for moder- 
ate k'^ three of the four intrinsic form factors resemble simple Gaussians, while the intrinsic 
neutron charge form factor requires an additional factor of k? in first approximation. Con- 
sequently, only a few terms of the Laguerre- Gaussian expansion are needed to obtain good 
fits, with higher-order terms used primarily for the estimation of the incompleteness error. 
Although it is possible to obtain fairly good fits using just 2 terms for Geu-, 4 for Ge-p-, or 6 
for Gwp and Gun-, in order to minimize model dependence and to evaluate incompleteness 
errors we employed 20 terms for each of the four form factors. 

The widths of the form factor bands are governed by the quality of the experimental 
data in the ranges of where data are available and for larger k'^ by the asymptotic 
envelopes indicated by dashed curves. Note that the uncertainties in the fitted form factors 
for k > kmax are reduced by the factor of 1/3 used in Eq. (|3l|) to transform from a uniform to 
a normal distribution and by the effect of constraints upon densities at large radii. Although 
the intrinsic form factors fitted to data for Gep, Gmp, and Gmu appear to decrease more 
rapidly than the k~^ envelopes, we prefer to employ these more generous uncertainties rather 
than to impose the steeper declines suggested by extrapolation from the measured into the 
unmeasurable region (where k > km)- The use of steeper envelopes would simply reduce 
the uncertainties in the extracted densities without affecting their central values. Therefore, 
k~^ envelopes were matched to fitted form factors at kmax = 6.5, 9.0, and 8.2 fm~^ for Gep, 
Gmp, and Gmu based upon the experimental Qmax for each form factor. However, uncritical 
application of the same procedure to Geu would suppress the high-frequency components of 
its intrinsic form factor too strongly because Qmax for Geti is presently too small to expect 
p{k) to decrease more rapidly than k~'^. Fig. |^ shows that the data presently available for 
Geu are compatible with the Galster parametrization, but the procedure used for the other 
form factors would cause Geh/Gd to decrease fairly rapidly beyond the range of these data. 
On the other hand, it is reasonable to expect Geh/Gd to decrease for Q"^ beyond a few 
(GeV/c)^, as observed in the other form factors. Therefore, in order to permit the positive 
slope for Geu/Gd to continue over a limited but larger range of Q^, we used the same value 
of kmax for both Geu and Gep even though the Geu data are limited to k < 5.4 fm~^. 
We believe that this compromise provides a more reliable extrapolation to higher and 
that the increased estimate of incompleteness error is more realistic, but obviously it is very 
important to acquire accurate data for Geu at higher Q"^ . 

These fits to intrinsic magnetic form factors do not change sign within the experimentally 
accessible region, k < km, but the fitted proton intrinsic charge form factor suggests a node 
between the present experimental limit, kmax = 6.5 fm~^, and ultimate limit, km = 9.5 
fm~^. Consequently, this model suggests a zero crossing in Gep near Q'^ ~ 10 (GeV/c)^. 
Figure | compares the form factor ratio /iGe/Gm deduced from the fitted form factors 
with experimental data using recoil polarization for the proton or using either recoil or 
target polarization for the neutron. For the proton we also show the linear parametrization 
proposed by Jones et al. for > 0.3 (GeV/c)^, while for the neutron we show a new 
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fit using the Galster parametrization, Eq. (|]), that gave A = 0.90 ± 0.02 and B = 3.8 ± 0.5. 
The data for the proton do not distinguish between linear and LGE parametrizations, but 
according to pQCD one would expect Gep/Gmp to approach a constant for sufficiently large 
Q^. Extrapolation of the LGE parametrization suggests that the asymptotic ratio will be 
very small, but data at much larger are needed to establish that level. An extension to 
9 (GeV/c)^ has been approved |7^, but larger remains desirable. Similarly, the present 
data for GEn/GMn are compatible with the Galster parametrization but remain limited 
to rather small Q^. Consequently, the extrapolation to larger is rather uncertain. If 
an approved experiment using the ^He(e, e'n) reaction |]73| achieves the proposed ±13% 
statistical uncertainty at = 3.4 (GeV/c)^, the error band will be reduced to about the 
same width and the extrapolation much improved. Nevertheless, there is little reason to 
expect the asymptotic limit to be reached earlier for the neutron than for the proton. 

Although a review of recent theoretical calculations is beyond the scope of the present 
work, it is probably worth mentioning a few which describe the new Gep/ Gup data relatively 
well. Among these the earliest is the chiral soliton model of Holzwarth pO| , which predicted 
the linear descrease with respect to and a sign change near 10 (GeV/c)2. More recently 



21| , modifications of the vector meson parameters were made to improve the fits to the 
neutron form factors, but the ratio Gmu/Gmp is not reproduced. Furthermore, because the 
chiral soliton model uses A^; = and = 1, Holzwarth found it necessary to artificially 
increase the soliton mass in order to obtain reasonable fits at large Q^. Alternatively, Lu et 
al. [0,^ obtained a good fit to the Gep/Gmp data for ^ 3 (GeV/c)^ by adjusting the bag 
radius in the cloudy bag model, but the ratio appears to level off well above the more recent 
data for higher Q^. Note that this model uses \e = Am = 1. The covariant calculation of 
Boffi et al. |76] using the point-form spectator approximation provides reasonably accurate 
predictions of the form factors for ^ 5 (GeV/c)^, although there remains a significant 
discrepancy for Gmp near the end of this range. The light-front calculations of Cardarelli 
and Simula using one-gluon exchange and the light-cone diquark model of Ma et al. |7S 
also reproduce the linear dependence of Gep/Gmp fairly well. 



B. Densities 

Proton charge and magnetization densities are compared in Fig. |^. Both densities are 
measured very precisely, with uncertainties at the origin better than 6% for magnetization 
or 8% for charge. Incompleteness dominates in the interior region while statistical errors 
become comparable in the surface region. As shown by the variation of Gep/Gmp in the top 
panel of Fig. ^, the new recoil-polarization data for Gep decrease more rapidly than either the 
dipole form factor or the magnetic form factor for > 1 (GeV/c)^. Consequently, we find 
that the charge density is significantly softer than the magnetization density of the proton. 
The densities obtained using LGE or FBE parametrizations are practically indistinguishable 
and are independent of the choice of h or Rmax over wide ranges. These densities are similar 
to the Gaussian densities one might expect from a quark model and are more realistic than 
the exponential density that results from naive nonrelativistic inversion of the dipole form 
factor. 

Neutron densities are shown in Fig. p. We find that the magnetization density for the 
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neutron is very similar to that for the proton, although the interior precision is not as good 
because the range of is smaller and the experimental uncertainties larger. Limitations in 
the range and quality of the Geu data presently available result in a substantially wider error 
band for the neutron charge density. Data at higher are needed to improve the interior 
precision, but a useful measurement of the interior charge density is obtained nonetheless. 
The positive interior density is balanced by a negative surface lobe. Note that polarization 
measurements are sensitive to the sign of the density. 

Whereas Figs. emphasize the interior densities, it is also of interest to compare these 
densities in the surface and tail regions. Figures use a factor of to emphasize these 
surface and tail densities. Although the densities are small, the reduced slopes seen between 
1 and 1.5 fm in the neutron magnetization and in both the charge and the magnetization 
densities for the proton are seen as significant peaks in r^p. Virtually identical features also 
emerge using the FBE parametrization. These features are independent of h for the LGE or 
Rmax for the FBE parametrization over wide ranges. Many attempts were made to suppress 
structures in r^p in the 1-1.5 fm region by limiting the number terms in the expansions 
or by application of stricter tail biases, but all modifications which did produce smoother 
r^p curves in this region damaged the fits to the form factor data for ^ 1 (GeV/c)^. 
Although it is difficult to prove that smoother fits do not exist, especially if one is willing to 
tolerate a moderate increase in x^, we were unable to produce acceptable fits without some 
structure in r^p in this region. On other hand, because local errors in momentum space 
can introduce artificial oscillations at large radii, we did apply an exponential tail bias for 
r > 2 fm where little structure is expected. Thus, the smaller oscillations for r > 1.5 fm 
are generally consistent with zero and can be suppressed using the tail bias with little effect 
upon the fits. We believe that the 2 fm matching radius is sufficiently large to have minimal 
infiuence upon densities for intermediate distances governed by data with of order several 
(GeV/c)2. 

The relatively small differences between Gun and Gu-p seen in Fig. ^ produce the small 
differences between neutron and proton magnetization densities shown in Fig. |. The peak 
of r^Pra is found at a slightly larger radius for the neutron than for the proton because the 
form factor decreases a little more rapidly with respect to . The secondary peaks in the 
1-1.5 fm region are also similar. Again, these comparisons are independent of the details of 
the analysis and are virtually identical using either LGE or FBE parametrizations. 

Therefore, the secondary r^p peaks in the 1-1.5 fm region appear to be essential features 
of the data rather than artifacts of analyses based upon linear expansions. While it is not 
possible to determine the physical mechanism for such features from data analysis alone, 
there is at least one simple candidate. The tensor interaction between quarks is expected 
to produce a small D-state component peaking at larger radius than the dominant S-state 
configuration, and the superposition of these components could yield a secondary peak at 
relatively large radius. 

C. Fitted Moments 

Moments and for each fit are listed in Table 0. The expansion coefficients require 
too much space to list but are available by request; also note that accurate reproduction of 



18 



the error bands would require full covariance matrices. Here we quote x^/^; where is the 
number of data points, because counting the number of degrees of freedom is not so clear 
when both high-A; and large-r constraints are applied. For each form factor we find that 
the 6 fits obtained for 3 possible choices of A used in either LGE or FEE parametrizations 
are essentially identical and give the same values of 'X^ /N. The normalization for Gep is 



consistent with unity within the 0.5% systematic uncertainty claimed by Simon et al. ||39 
for their data at low Q^. For the Gmp at low we employed the results of Hohler et al. 
who adjusted the relative normalizations of several data sets to a common standard. The 
normalization produced by the present fit is consistent with the systematic uncertainty in 
that standard. Except for Gmti, data selected from several sources appear to be mutually 
consistent and the quality of the fitted form factors is very good. Although the low-Q^ data 
for Gmu have improved in recent years, significant systematic discrepancies remain. Recent 
data from Refs. |53| , ^ , ^ , |50| with small statistical uncertainties suggest a small dip near 0.2 
and a peak near 1 (GeV/c)^. However, the data from Refs. PU| , ^ are inconsistent with the 
fit and inflate x^- Nor do the data for Gmu reach sufficiently low to strongly constrain the 
normalization and the data sets are not entirely consistent either; consequently, the present 
analysis suggests that a 2% normalization error remains. We decided to retain the data from 
Refs. pB|^ , despite their deviation from the fit because we are not entirely convinced of 
the transportability of the efficiency calibration from a hadron facility to an electron facility. 

The final two columns of Table |I| list the Sachs radius, C,, defined by Eq. (|2T| ) in terms 
of the initial logarithmic derivative of the Sachs form factor and the transition radius, R, 
obtained from integral moments of the density fitted for a specific A according to Eq. (p5|). 
The dependencies of rms radii for Gep, Gmp, and Gmu upon A are consistent with Eq. (p2D, 
showing a significant discrete ambiguity arising from the model dependence of the form 
factor to density inversion. By contrast, for each of these form factors all 6 determinations 
of C, are consistent with each other, demonstrating that the fitted is a model-independent 
property of the Sachs form factor. 



1. Proton charge radius 



The analysis of Simon et al. has been accepted for about two decades as the definitive 
determination of the proton charge radius, but more recently there has been renewed interest 
in that quantity now that the finite-size corrections have become the dominant uncertainty 
in theoretical calculations of the IS Lamb shift in hydrogen. For example, Melnikov and 



van Ritbergen |^ argue that the Lamb shift provides the most accurate measurement of 
the proton charge radius and deduced a value 0.883(14) fm that is somewhat larger than 
the 0.862(12) fm obtained by Simon et al. Although our definition for the intrinsic charge 
radius depends upon the choice of A^; employed to fit the Sachs form factor, the definition 
generally employed by QED theorists corresponds to the quantity we labelled as ^Ep that 
is based upon the initial slope of GEpiQ"^) and is independent of A^;. Thus, it is the fitted 
value of ^Ep that should be compared with the Lamb shift result. 

Our fit using the LGE parametrization is compared with the data from Simon et al. in 
Fig. ^. Also shown is the monopole fit made by Simon et al. to the data for < 2.3 (GeV/c)^ 
that were available at that time. That analysis gave a smaller value, ^Ep = 0.862 ± 0.012 
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fm, but does not fit the low data as well as our LGE fit. Our fit employs the entire data 
set described in Sec. Ill Q even though Fig. ^ shows only the lowest region. Fits made 
using different values of A^; or using the FBE parametrization are indistinguishable and give 
values for that are consistent within their quoted uncertainties. Therefore, we claim 
that = 0.88 ±0.01 fm represents a model-independent property of the experimental data 
even if its interpretation as a charge radius depends upon the choice of A^- This value is 
consistent with the rms radius derived by Melnikov and van Ritbergen |]7D[ from the Lamb 
shift. 

Coulomb distortion may also affect the charge radius obtained by electron scattering. 
Rosenfelder analyzed this effect using a distorted wave Born calculation to correct mea- 
sured electron scattering cross sections for Coulomb distortion, thereby obtaining effectively 
plane-wave cross sections. The charge radius was then obtained by fitting the resulting form 
factors with low-order polynomials. Although the form factor corrections were typically less 
than 1%, the value of ^Ep extracted from the adjusted data of Simon et al. nevertheless 
increases by about 0.008 - 0.013 fm depending upon the fitting strategy and the degree 
of the polynomial. However, that analysis employed a Coulomb potential obtained from a 
charge density of the form 

This form not only lacks the Lorentz-contraction factor for spatial frequency, but uses a 
value for A^; = —1/2 that is inconsistent with the high behavior of the form factor. The 
question of Coulomb distortion of the low form factors merits further investigation, but 
a consistent relativistic relationship between form factor and density must be employed. 
Nevertheless, the magnitude of that correction appears to be smaller than the uncertainty 
in the fitted quantity. 



2. Proton magnetization radius 

With S,Mp = 0.85 ± 0.03 fm the proton magnetization radius appears to be slightly 
smaller than its charge radius, as expected from Fig. ^ but the uncertainty is as large as 
the difference because the data at very low Q"^ are less precise for Gmp than for Gep- A 
substantial part of the uncertainty in ^Mp is due to the uncertainty in normalization. If we 
constrain Mq to unity, the changes in fitted form factors and densities are relatively small, 
but we obtain values for ^Ep = 0.862 ± 0.006 fm and ^Mp = 0.835 ± 0.006 fm that appear 
to be much more precise. The constrained fit to Gep is close to the result of Simon et al. 
shown in Fig. ^ at very low Q^. This analysis demonstrates that there is an appreciable 
difference between the proton charge and magnetization densities, but it also highlights the 
importance of precise absolute normalization at very low Q^. 



3. Neutron magnetization radius 



Our fits to the Gmu data give a value for ^Mn = 0.92 ± 0.07 fm that has substantial un- 
certainty because the lack of data for very low permits significant widening of the error 
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band as 0. By contrast, Kubon et al. [^] obtained a value ^Mn = 0.873 ± 0.011 fm 

that appears to be much more precise. However, their continued-fraction parametrization 
automatically constrains the normalization at = 0. If we constrain the normalization 
by requiring Mq = 1, then the LGE analysis gives a value for C,Mn = 0.881 ± 0.018 fm that 
agrees with Kubon et al. Note that we included data from Refs. [0,0 that were omitted 
by Kubon et al. and that deviate strongly from our LGE fit, but these data appear to have 
little influence upon the fitted normalization and radius. Although the fitted normalization 
is consistent with unity, the 3% uncertainty does have an appreciable effect upon the uncer- 
tainty in ^Mn- We chose not to constrain the normalization in the standard analysis because 
the systematic errors in neutron efficiency have been a big problem for Gmu measurements 
and there remains significant scatter among recent experiments at low Q^. 



4- Neutron charge radius 



The charge radius for the neutron can be expressed in terms of Dirac and Pauli form 



factors as 



-6 



3k„ 
2m^ 



(38) 



where the first term is sometimes described as the intrinsic radius {e.g. [^,^) while the 
second term is called the Foldy term and is attributed to a charge separation induced by the 
zitterbewegung motion of the magnetization density. We prefer to describe the first term as 
the Dirac radius because it is derived from the initial dependence of the Dirac form factor 
and to reserve the term intrinsic radius for Eq. (]TP|), which is based upon a moment of the 
radial function that we identified as the intrinsic charge density. However, for the present 
purposes it will be clearer to refer to the radius based upon the Sachs form factor as the Sachs 
radius. The observation that the Foldy term, equal to -0.126 fm^, is by itself almost equal 
to the mean-square charge radius obtained from bm has generated considerable discussion. 
Furthermore, as discussed in Sec. |111 C| , a substantial disagreement remains between the 
results from Dubna and those from Oak Ridge, Garching, and Argonne. For example, 



Alexandrov |^6[ argues that the Foldy term should be discarded and that the pionic cloud 
should make the mean-square Dirac radius negative. He further claims that experiments 
giving bne ~ —1.31 fm are likely to suffer from serious experimental or interpretative errors 
because the corresponding mean-square Dirac radius would be positive. On the other hand, 
Glozman and Riska |^| calculated that the pion loop contribution to the Dirac radius is 
negligible. Using a Foldy- Wouthuysen analysis of the interaction of a particle with internal 
structure with an external electromagnetic field, Bawin and Coon |^ demonstrated that 
the Foldy term is cancelled by a higher-order term arising from the Dirac form factor, leaving 
the Sachs radius as the dominant coefficient, independent of the dynamics responsible for 
the neutron form factors. 

From a more microscopic point of view, Isgur used a quark model to argue that the 
observation of a very small Dirac radius is a potentially misleading accident and that the 
relativistic boost cancels the Foldy term such that the slope of Geu does provide the sec- 
ond moment of the intrinsic charge distribution [Q. Cardarelli and Simula [§3|,|7^ showed 
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that this cancellation depends upon neglect of transverse momenta and that quark spin-spin 
interactions that break SU(6) symmetry provide a mixed-symmetry S' component that en- 
hances Geti and provides a fairly accurate fit to the recoil-polarization data for fJ^pGEp/GMp- 
Leinweber et al. argued that chiral perturbation theory provides model-independent 
constraints on the dependencies of nucleon magnetic moments and charge radii upon quark 
masses that demonstrate that the similiarity between the Sachs radius and the Foldy term 
is purely accidental. 

Our fitted mean-square radius for the neutron charge density is largely determined by 
and is completely consistent with the datum of Ref. that was included in the analysis. 
This quantity is free of discrete ambiguity because Mq vanishes and the analysis procedure 
enforced that constraint explicitly. Unfortunately, the electron-scattering data are not suf- 
ficiently precise at very low to resolve the controversy concerning the sign of the Dirac 
radius. In the absence of a datum for M2, the fit to the data for Geu gives M2 = — 0.187±0.04 
fm, which is about one or two standard deviations larger in absolute magnitude than the 
Dubna or Oak Ridge results, respectively, but is much less precise than either — the error 
bands on Gsn with or without the M2 datum overlap almost completely, even at low Q^. The 
effect of the constraint on M2 is shown in Fig. |10|, where the shaded band for Fin represents 
the unconstrained fit to Geu data that has a negative Dirac radius while the cross-hatched 
band with a positive Dirac radius includes the M2 datum and is consequently narrower as 
0. The rather small difference between these fits is confined to ^ 0.2 (GeV/c)^, 
but because Fin < over most of its measured range the fit with negative Dirac radius 
requires a sign change near ^ 0.14 (GeV/c)^. Therefore, the present electron-scattering 
data offer very little sensitivity to M2. Furthermore, because the nuclear physics corrections 
needed to extract Geu for low are substantial even for polarization methods, it is unclear 
whether one can ever expect better accuracy from nuclear physics than atomic physics mea- 
surements of M2. We consider neither the theoretical argument for negative Dirac radius 
nor the limited experimental evidence for a sign change in Fin at low compelling. 



V. DISCUSSION 

A. Comparison with Nonrelativistic Analyses 

Several analyses have appeared recently in which the parameters of the Galster model 
were fitted to selected data for Geu at low and a density extracted using the nonrela- 
tivistic inversion formula given by Eq. (^. Examples of this type are shown in Figs. pJ] - |T^ . 
Figure O compares fits of this type to recent data using the d(e, e'n), d{e, e'n), or ^He(e, e'n) 



reactions. The slope of the form factor obtained from the Oak Ridge value for bm is shown 



as a line segment. The result obtained by Platchkov et al. [q^ from an analysis of elastic 
cross sections for electron scattering from deuterium using the Paris potential is shown as 
the dashed curve and lies well below the data obtained from polarization measurements. 
However, variations of ±50% in Geu were found using different realistic nucleon- nucleon 
potentials; the result using the Paris potential is quoted most often, but the result using 
the Argonne V14 potential is closer to the modern data. The dash-dotted curve shows a fit 
by Schmieden |B6| to the data from Mainz, excluding the point for ^He(e, e'n) at = 0.4 
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(GeV/c)^; this fit gives the highest values for Geh- The original Galster model is shown as 
a solid curve and our fit based upon Eq. (|) is shown as a dotted curve and lies between the 
Paris and Mainz results. The fit by Zhu et al. [0] is also close to the dotted curve. Note 
that our fit uses the entire data set described in Sec. [Ill Q , including the slope at the origin, 
while the Mainz fit used a smaller subset. Even though the data employed by Galster et al. 
had much larger uncertainties, their result is remarkably close to the present analysis. 
However, the apparent agreement of the original Galster parametrization with more modern 
data must be judged as fortuitous. 



Buchmann et al. [5B,5U| argue that the neutron charge density refiects differences between 
the spatial distributions of constituent quarks induced by the color hyperfine interaction and 
that because the same interaction is responsible for the N — A mass splitting the N ^ A 
quadrupole form factor is related in a simple manner to Geu- In Ref. they used the 
Galster parametrization to fit a selection of the Geu data and used nonrelativistic inversion 
to obtain a density similar to the solid curve in Fig. |T^. 

The nonrelativistic densities obtained from these fits are compared in Fig. to the 
present relativistic LGE result. Here we chose A^; = to minimize the differences between 
relativisitic and nonrelativistic inversion formulas; however, it is important to remember that 
the Galster form factor is inconsistent with the relativistic inversion procedure for A^; = 0. 
The upper panel, emphasizing the interior density, demonstrates that the naive Fourier 
transform tends to produce a rather hard core and an unphysical cusp at the origin. The 
relativistic transformation, by contrast, softens the interior density and eliminates its cusp, 
providing a much more plausible charge density. The lower panel, emphasizing the surface 
lobe, shows that all three nonrelativistic densities have positive peaks at smaller radii and 
less surface charge than our result. 



B. Approach to scaling 



The asymptotic behavior of the fitted form factors is illustrated in Fig. 13 by mul- 
tiplication by Q^. The uncertainties in the fitted form factors are clearly dominated by 
experimental uncertainties where data are available, while the expansion of the error bands 
for larger is governed by the incompleteness errors in the range k^ax < k < km- The 
scaling behavior of Gup appears fully developed because the data up to Ql^ax = 31 (GeV/c)^ 
reach a sufficiently large ratio kmax/km = 0.95 to strongly constrain the asymptotic limit. 
Although the fit to Gmti is compatible with pQCD scaling, the uncertainties at large grow 
much more rapidly because here kmax/km = 0.86 is smaller and the data for kmax/km > 0.79 
are much less precise. Furthermore, because scaling is not fully developed for Gmp until 

^ 20 (GeV/c)^ we should expect that data for Gmu at higher will be needed to es- 
tablish its asymptotic limit. The data for Gep do not show scaling behavior for Qmax = 5.54 
(GeV/c)^, where kmax/km = 0.78 — despite relatively large uncertainties in extrapolation, 
the present fit suggests a sign change near ~ 10 (GeV/c)^. The new recoil polarization 
data suggest that Gep and Gmp differ dramatically for ^ 3 (GeV/c)^ and it is clearly 
crucial to extend the Gep data to higher Q^. It will also be important to check those data 
using another technique. The results of an improved Rosenbluth experiment are expected 
soon 143 
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The present data for Geu are too limited in both range and precision to address the 
question of scahng. High-precision data for < 1.5 (GeV/c)^ from a recent d(e, e'n) 
experiment are expected soon and an approved proposal for ^He(e, e'n) should extend 



the range to 3.4 (GeV/c)^ in a couple of years [^]. However, experience with the other three 
form factors suggests that one must approach 20 (GeV/ c)^ to determine the asymptotic limit. 

Quark helicity conservation suggests that Q^F2/Fi should approach a constant in the 
asymptotic limit of large Q^. The Dirac and Pauli form factors are related to Sachs form 
factors by 

F^ = ^^±1^ (39a) 



1 + r 

Gm — Ge /Qnu^ 
kF2 = (39b) 



and their ratio is given by 



Fi T + g 



(40) 



where g = Ge/Gm can be measured directly by either recoil or target polarization. Recoil- 
polarization data for this ratio are compared in Fig. ^ with bands constructed from the 
present fits to Sachs form factors. Gayou et al. [|l6l observed that the proton recoil polar- 
ization data appear to reach a plateau in the range 2 < < 6 (GeV/c)^ when scaled by Q 
instead of the expected . Ralston et al. [0,^ suggested that orbital angular momentum 



in the quark distribution could explain an asymptotic behavior of the form F2/F1 oc Q ^. 



Later Miller and Frank argued that substantial violation of quark helicity conservation 
should be expected for intermediate when Poincare invariance is imposed upon relativistic 
constituent quark models. 

Recognizing that g is small compared with r for > 6 (GeV/c)^ and that the model 
imposes constraints upon the Sachs form factors for high Q'^ that inhibit the growth of 
the error band for F2/F1, we can extrapolate F2/F1 beyond the present experimental range 
for Gep- This extrapolation, shown in Fig. |15|, shows that the data are consistent with 
quark helicity conservation for ^ 20 (GeV/c)^. Therefore, although the present model is 
consistent with the observation by Gayou et al. that QF2/F1 is approximately constant in 
the range 2 < < 6 (GeV/c)^, we attribute that observation to the appearance of a broad 
maximum in QF2/F1 rather than to the onset of true scaling. 

C. Quark densities 

Assuming isospin symmetry and neglecting strange quarks, nucleon charge densities can 
be expressed in a simple two-flavor quark model as 

4 1 

Ppir) = -u{r) - -d{r) (41a) 
2 2 

Pnir) = --u{r) + -d{r) (41b) 
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where u{r) is the radial distribution for an up quark in the proton or a down quark in the 
neutron while d{r) is the distribution for a down quark in the proton or an up quark in the 
neutron. Thus, the quark densities are obtained from nucleon charge densities using 

u{r) = pp{r) + ^p„(r) (42a) 
d{r) = pp(r) + 2p„(r) (42b) 

where u{r) and d{r) are normalized to unity according to 

POO 

/ dr r'^q{r) = 1 (43) 
Jo 

but need not be positive everywhere. There is no guarantee that these combinations of radial 
densities obtained from form factor data by relativistic inversion must be positive, nor are 
the densities derived from positive-definite matrix elements. Within the quark model one 
could decompose the densities 

qir) = q^{r) + qs{r) - qs{r) (44) 

for each flavor, q, into valence (v) and sea (s) contributions which would be expected to 
be positive definite but cannot be separated using only Sachs form factors. The valence 
distributions, q^, are normalized to unity while the sea distributions, and g^, must have 
equal normalizations. Thus, within this picture the contributions from quarks should be 
positive but antiquarks in the sea could produce regions where u{r) or d{r) might become 
negative. 

The quark densities obtained with LGE densities for A^; = 2 that are displayed in Fig. 
p!6| are determined with relatively small uncertainties and are predominantly positive, as 
expected. We find that the pair of like quarks has a somewhat broader distribution than 
that of the unlike quark, being depleted in the interior and enhanced at the surface. In 
this model, the neutron charge density arises from incomplete cancellation between charge 
densities for up and down quarks, resulting in positive core and negative surface charges. 
The broader distribution for like quarks is consistent with the repulsive color hyperfine 
interaction between like quarks needed to explain the — A mass difference. This picture 
is also consistent with the model of a pion cloud surrounding a three-quark core, but in 
that model one might expect to find a slightly negative d{r) near the surface due to the 
antiquark content of the pion. The present data are marginally consistent with a slightly 
negative d{r) near 1.0 fm, but more accurate data for the neutron charge density would be 
needed to reduce the uncertainty in d{r) before drawing a definitive conclusion. 



D. Discrete Ambiguities 

Although discrete ambiguities in A^; and Am do not affect fitted form factors in the 
range where data are available, the choice of A does affect the growth of the uncertainties 
in extrapolated form factors as increases beyond the measured range. For example. Fig. 

shows form factors fitted using the LGE parametrization with A^ = and Aa/ = 1 as 
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suggested by the relativistic soliton model. Although the relative uncertainties become quite 
large for beyond the range of the experimental data, the uncertainties in the form factors 
actually remain small because, with the exception of Geu-, the form factors are greatly 
reduced as their experimental limits are reached. Consequently, the contribution of the 
uncertainties in form factors at large upon uncertainties in radial densities is relatively 
insensitive to the choice of A. However, the change in fitted density due to a change of A need 
not be contained within the fitted error band — the bands do not accommodate discrete 
ambiguities in the model. 

The sensitivity of fitted charge densities to the choice of A^; is illustrated in Fig. ff^. These 



figures were made with the LGE parametrization, but very similar results are obtained with 
the FEE parameterization. The smoothest results at large radii are obtained with A = 0, 
whereas larger values of A tend to pull the density inward and to amplify oscillations at large 
radii. This behavior can be understood by interpreting Eq. ([l^) in terms of the convolution 
theorem for Fourier transforms. Expressing this equation in terms of 

m = {I - ^)-'G{-^) (45) 

one finds that p{k) for A > is obtained from the Lorentz contracted form factor by decon- 
volution of a resolution function with mean square radius equal to 3A/2m^. This resolution 
function originates in the zitterbewegung and is characterized by the nucleon Compton wave- 
length. As discussed in Sec. [IV B| , acceptable fits to the Sachs form factor data for ~ 1 



(GeV/c) using A = 2 seem to require structure in the radial densities in the 1-1.5 fm region. 
Reducing A tends to smooth out such structures, but sacrifices the high limit. Therefore, 
although A = 2 provides the most natural implementation of pQCD scaling, accurate repro- 
duction of the data using that representation of the intrinsic form factor appears to require 
oscillations in the radial density with a wavelength of order 0.7 fm. 

Despite this ambiguity in the relationship between form factor and density, the qualitative 
conclusion that u{r) is broader than d{r) depends only the assumption of isospin symmetry 
and the observation that the neutron charge density is positive in the interior and negative at 



its surface. The quark distributions derived using A^; = shown in Fig. |T9| are qualitatively 
similar to those shown for A^; = 2 in Fig. but are slightly more diffuse. The choice A^; = 
is appropriate for the soliton model, is consistent with nonrelativistic expectations for small 
Q^, and is favored by the radius obtained from the Lamb shift. However, it appears to be 
inconsistent with pQCD. Therefore, in the absence of a unique relationship between form 
factors and densities, it appears necessary to select the appropriate value of A based upon 
the intended application. For long- wavelength properties one should use A^ = 0, but for 
extrapolation to the pQCD limit one should employ A^; = Am = 2. 



E. Alternative parametrizations 



Without a unique relationship between form factor and intrinsic density, one may ques- 
tion the value of densities extracted by the present techniques. It is clear that simple 
parametrizations like the dipole model cannot fit the proton or Gmu data over a wide range 
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of Q , nor is it likely that the Galster parametrization will continue to fit Geh at higher Q . 



The empirical parametrization proposed by Bosted ^4 



G oc (1 + aiQ + aaQ' + a^Q^ + a^Q^^ (46) 

fits the data for large well and is consistent with pQCD, but its odd powers of Q are 
incompatible with the interpretation of the form factor as the Fourier transform of a radial 
density and with the moment expansion for small Q^. Furthermore, it is not sufficiently 
fiexible to provide realistic error bands, especially if the odd powers are eliminated. By 
inclusion of a term. Brash et al. also sacrificed the pQCD limit in order to improve 
the quality of the fit for finite Q^. 



Kubon et al. fit a subset of the Gmh data using a continued-fraction parametrization 
of the form 

GMnm = — ^ (47) 

carried to fifth order. This parametrization provides a good fit to the data for < 4 
(GeV/c)^ using 5 parameters, but the parameters do not decrease with order and the fit 
depends upon fairly delicate cancellations. Adding additional terms to extend the range of 
changes the lower terms. Furthermore, this fit does not conform with the asymptotic 
Q""^ behavior expected by pQCD unless a fairly complicated constraint of the form 

M5 + &2(&4 + &5) =0 (48) 

is imposed to eliminate the contribution. The constraint increases in complexity as 
additional terms are included. A comparison between the Kubon parametrization and our 
LGE fit with Xm = 2 is shown in Fig. EUl The Kubon analysis included only the data 



indicated by filled circles and was limited to < 4 (GeV/c)^ while our analysis used 
all data shown and extended to 10 (GeV/c)^. The LGE parametrization fits the data 
well over a broader range and is compatible with pQCD, whereas the continued- fraction 
parametrization behaves badly soon after the range fitted by Kubon et al. By not imposing 
the pQCD constraint, the extrapolation deteriorates quite quickly. We believe that the 
continued-fraction method also underestimates the uncertainty in the rms radius due to the 
strong correlations among its parameters and its built-in normalization constraint at = 0. 

A rather different phenomenological parametrization can be made in the context of 
the vector meson dominance model at modest matched to pQCD at large Q^, denoted 
VMD-|-pQCD. That approach was pioneered by Gari and Kriimpelmann |p5| , |96| , |97| , |98[| and 
recently refined by Lomon [|99| , |100|| . Similarly, the classic dispersion-theory analysis of Hohler 



et al. |2^ has recently been updated by Mergell et al. ||2^ to handle better the requirements of 
unitarity and the approach to the pQCD limit. These approaches have the advantage that all 
four electromagnetic form factors are analyzed simultaneously, thereby relating their isospin 



structure to an underlying model, and can be extended to timelike momentum transfer p6 
By contrast, our approach is limited to spacelike momentum transfer and must construct 
the isospin form factors from four independent fits to individual form factors. Both the 
VMD and dispersion theory approaches appear to be capable of fitting the data as well 
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as our linear expansion analysis, although the data have improved considerably since the 
analysis of Ref. |25]. However, we omit detailed comparisons here because these models do 
not consider radial densities. 



F. Importance of Geu data at higher 

The present data for Geu do not extend high enough in to determine the interior 
charge density as accurately for the neutron as for the proton or to permit reliable extrap- 
olation to the scaling regime, but new data expected from an approved proposal at 



Jefferson Laboratory should help considerably. The impact of extending the Q range to 
3.4 (GeV/c)^ is illustrated in Fig. ^ This analysis was performed using A^; = 0, which 
permits the greatest latitude at high Q^. The left column shows the form factor and den- 
sity fitted to published data (shown by open symbols), while the middle and right columns 
show the effect of pseudodata (shown as filled symbols) for two hypothetical scenarios. The 
middle scenario assumes that the new data would follow the Galster parametrization while 
the right scenario assumes that the new data would fall more rapidly than the dipole form 
factor for > 2 (GeV/c)^. Both scenarios are compatible with the uncertainties extrapo- 
lated from the fit to the present data and can be fitted well, with error bands for ^ 4 
(GeV/c)^ consistent with the anticipated experimental precision. The reduction of the in- 
completeness error obtained by extending the measurements to higher greatly improves 
the precision of the interior charge density. If measurements at higher come close to the 
Galster parametrization, then the error band would be reduced in width with little change 
in its centroid. On the other hand, if new measurements of GEn/Go decrease with in a 
manner similar to the proton charge form factor, the softer charge density would be reduced 
in the interior, moving toward the lower edge of the present error band. These scenarios 
have quite different asymptotic values for Q^Geti, but the present data cannot distinguish 
between them. Furthermore, the recent VMD+pQCD analysis by Lomon [|100|| suggests that 



Geti/Ge could reach an asymptotic value substantially higher than predicted by the Galster 
parametrization. Therefore, it is very important to extend Geu data as far as possible in 



VI. SUMMARY AND CONCLUSIONS 

We have employed expansions in complete sets of radial basis functions to parametrize 
nucleon Sachs form factors in terms of charge and magnetization densities. Our selection 
of data emphasizes recent polarization data. The inversion from form factor to density is 
based upon relativistic models in which the spatial frequency k = Q/y/l + t in the rest 
frame is related to the momentum transfer Q in the Breit frame by Lorentz contraction. 
The maximum possible frequency sampled by electron scattering is then k < where 
km = 2m is determined by the nucleon Compton wavelength. A variety of models produce 
inversion formulas of the form 
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but differ in the choice of A. By considering the asymptotic form of G{Q'^), we can Umit 
the exponent A to 0,1,2. The relativistic soliton model suggests {A^; = O.Xm = l}i the 
original quark cluster model suggested Xe = Am = 1, and a more symmetric version of the 
quark cluster model gives choices A^ = Am = 2 that are compatible with pQCD without 
the somewhat artificial constraints upon p{km) needed by the other models. In most of 
this paper we parametrized the Sachs form factors using the Laguerre-Gaussian expansion 
(LGE) and derived densities using Ag = Am = 2, but we have also analyzed the impact of 
the discrete ambiguity in A upon the radial densities. Although some of the details of the 
radial densities are affected by the discrete ambiguity in the relativistic inversion formula, 
their qualitative features are independent of A. 

We find that virtually identical fits to the Sachs form factors are obtained with either LGE 
or Fourier-Bessel expansions (FEE) and within a wide range these fits are independent of 
details of the parametrization, such as number of terms, radial scale parameters, and tail bias. 
The fitting procedure uses p{k) pseudodata for k > k^ax to estimate the incompleteness error 
in radial density due to the limitation of experimental data for the range k < k^ax < km- 
For a given choice of A, the radial densities fitted with either LGE or FEE expansions are 
practically identical and are insensitive to details of the analysis. We find that the proton 
charge density is significantly broader than its magnetization density, consistent with the 
observation in recent recoil-polarization measurements that Gep/Gmp decreases in an almost 
linear fashion for 1 < < 6 (GeV/c)^. Our result for the proton Sachs charge radius is 
consistent with a recent determination based upon the IS Lamb shift in hydrogen. Similarly, 
we find the magnetization density is slightly broader for the neutron than for the proton. 
Each of these three densities exhibits a secondary peak in r^p(r) near 1 — 1.5 fm that 
cannot be suppressed without sacrificing the fit to the corresponding form factor for > 1 
(GeV/c)^. This structure appears somewhat stronger for A = 2 than for A = 0, but A = 2 
provides the clearest extrapolation to the pQGD limit. 

The recent recoil-polarization data for F2p/Fip appear to favor scaling with rather 
than the Q^^ expected from quark helicity conservation in pQGD. Although that observation 
has stimulated some speculation about violations of quark helicity conservation due to orbital 
angular momentum or imposition of Poincare invariance, in our analysis with Xe — Am = 2 
we find that although QF2p/Fip appears to be nearly constant for 2 < < 6 (GeV/c)^ 
we nevertheless obtain a constant asymptotic value for Q'^F2p/Fip for ^ 20 (GeV/c)^ 
where Gmp scales with Q^"^. Therefore, we find that the data are consistent with a broad 
maximum in QF2p/Fip and do not require true scaling. 

We have compared the LGE parametrization for Geh to fits based upon the Galster 
parametrization. Although the range of remains too small to discriminate between these 
models, the Galster parametrization cannot be inverted using a relativistic relationship 
between intrinsic spatial frequency and Breit-frame momentum transfer unless A > 1. The 
traditional nonrelativistic inversion of the Galster form factor produces a charge density with 
an unphysical cusp at the origin while the relativistic fit to the data using the LGE form 
factor softens the interior density and removes the cusp. However, the incompleteness error 
in the neutron charge density remains fairly large because the available Geh data are limited 
to small and the data for 0.4 < < 1.6 (GeV/c)^ have relatively large uncertainties. 
More precise data for 0.5 < < 1.5 (GeV/c)^ are expected soon and an experiment for 
Q'^ < 3.4 (GeV/c)^ is in preparation. These data should improve the accuracy of the neutron 
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charge density considerably, but data approaching 20 (GeV/c)^ will probably be needed to 
test scaling in the neutron and in the isospin form factors. 

Combining the neutron and proton charge densities, we deduced the up and down 
quark radial distributions assuming isospin symmetry and neglecting heavier quarks. This 
schematic model suggests that the distribution is slightly broader for up quarks than for 
down quarks in the proton. With A^; = 2 we also observe a statistically significant negative 
density for down quarks near 1 fm that might be attributed to the d content of the pion 
cloud. 

Although we cannot claim that there is a unique relationship between form factors and 
densities, expansion of densities in a complete radial bases provides physically appealing 
parametrizations of form factor data that are applicable over a wide range of Q^. The use of 
linear expansions in complete bases minimizes the model dependence of the fitted form fac- 
tors and provides more realistic error bands in both spatial and momentum representations. 
Therefore, even if the identification with static densities is discounted, the fitted densities 
do provide useful parametrizations of the form factors nonetheless. The choice = Am = 2 
automatically satisfies pQCD scaling and provides a natural means for extrapolating form 
factors to higher Q'^ for the purpose of planning future experiments. 
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FIGURES 




FIG. 1. The ratio between Sachs form factor with A = 2 and the dipolc form factor is shown 
for a Gaussian intrinsic density using several values of the oscillator parameter, b, listed with units 
in fm. 




[(GeV/c)2] Q2 [(GeV/c)2] 

FIG. 2. The bands show fits to selected data for nucleon electromagnetic form factors using 
the LGE parametrization with A^; = Am = 2. For Geu the solid line shows a two-parameter fit 
based upon the Galster parametrization. 
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[(GeV/c)2] [(GeV/c)2] 

FIG. 3. Fourier transforms of the nucleon charge and magnetization densities are shown as 
error bands. The LGE parametrization was used with = Am = 2. The dashed hues show the 
upper Hmits used for estimation of incompleteness errors. The vertical lines divide the regimes of 
spacelike and timelike Q^, where the timelike threshold, < —Am?, is approached in the limit 
k"^ oo. Note that electron scattering is limited to the spacelike regime, wherein — >■ oo is 
represented by A; — ^ 2m. 
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FIG. 4. Comparison between data for Ge/Gm obtained from polarization measurements with 
fits made to the entire data sets employed for nucleon electromagnetic form factors. Results for the 
LGE parametrization A^; = Am = 2 are shown as bands. Also shown are the linear parametrization 



proposed by |14] for the proton and a fit based upon the Galster parametrization for the neutron. 
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FIG. 5. Comparison between charge {pch) and magnetization {pm) densities for the proton 
fitted using the LGE parametrization with A^; = Am = 2. Both densities are normalized to 
/ dr r'^p{r) = 1. 
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FIG. 6. Charge (pch) and magnetization (pm) densities for the neutron fitted using the LGE 
parametrization with Xe = Am = 2. 
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FIG. 7. Comparison between proton charge and magnetization densities using a factor of to 
emphasize the surface and tail regions. The fits used the LGE parametrization with Xe = Xm = 2. 
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FIG. 8. Comparison between neutron and proton magnetization densities using a factor of 
to emphasize the surface region. The fits used the LGE parametrization with Am = 2. 
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FIG. 9. Low fits to proton charge radius. The solid hne shows our LGE fit with A^; = to 



the entire data set while the dashed line shows the monopole parametrization of Simon et al. [39|. 
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FIG. 10. Sensitivity to the Dirac radius: the hatched (dotted) bands were obtained from fits 
to the Geu data that include (omit) the Oak Ridge datum for the neutron charge radius. 
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FIG. 11. Selected data for Geu at low are compared with fits based upon the Galster model. 
The solid curve is the original Galster fit while the dotted line is a new fit based upon the entire 
data set considered in this work. The dashed curve is a fit by Platchkov et al. [65| to elastic 
scattering by deuterium based upon the Paris potential. The dash-dotted curve was fitted to a 
subset of the Mainz data by Schmieden p6|. 
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FIG. 12. Neutron charge densities obtained from nonrelativistic Fourier transform of fits using 
the Galster model are compared with the present results (band) using the relativistic transformation 
with Xe = 0- 





1 2 
[(GeV/c)^] 



5 10 
[(GeV/c)2] 



40 



FIG. 13. The approach to scahng is shown by multiplying Sachs form factors by Q^. For Geu 
the original Galster model is shown as a dashed curve and a new fit as the solid curve. The LGE 
parametrization was used with A^; = Am = 2. 
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FIG. 14. Data for F^/Fi are compared with the usual scaling expected from pQCD or with 
Q^^ scaling recently proposed by several authors. Fitted bands employ the LGE parametrization 



with A^; = Aa/ = 2. The data are shown as squares |101], triangles circles [14,16|. 
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FIG. 15. Data for F2/F\ are compared with the usual scaHng expected from pQCD or with 
Q^^ scahng recently proposed by several authors. Fitted bands employ the LGE parametrization 
with \e = Am = 2. 
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FIG. 16. Quark densities obtained from proton and neutron charge densities using the LGE 
parametrization with Ae = 2. 
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FIG. 17. The bands show fits to selected data for nucleon electromagnetic form factors using 
the LGE parametrization with Ag = and Am = 1 as suggested by the soliton model. 
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FIG. 18. Discrete ambiguities in the charge densities. The value of A^,' increases from to 2 
from left to right. The curves for A^; = are reproduced, without error bands, in the middle and 
right columns for comparison. 
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FIG. 19. Quark densities obtained from proton and neutron charge densities using the LGE 
parametrization with A^ = 0. 
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FIG. 20. Comparison with the Gmu analysis by Kubon et al. ||5^. The band shows our LGE fit 
with Am = 2 to the entire data set while the solid line shows the continued-fraction fit by Kubon 
et al. (5^ to a subset of the data. The data included (omitted) by Kubon et al. are indicated by 
filled (open) circles. 
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FIG. 21. The sensitivity of the neutron charge density to extension of the experimental range 
of is ilhistrated by comparing a fit to published data (open circles) with two scenarios that 
include pseudodata (filled circles) for = 2.4,3.4 (GeV/c)^. The middle column assumes that 
new data would follow the Galster parametrization, shown by the solid curve, while the right 
column assumes that Geu/Gd would decrease for > 2 (GeV/c)^. Densities for the Galster 
model use nonrelativistic inversion. All fits use the LGE parametrization and relativistic inversion 
with A^; = 0. 
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TABLES 



TABLE I. Data selection summary. 



quantity reaction (GeV/c)^ Ref. 



method 



Gmv P{e, e') 



0.017 - 0.15 
0.16 - 31.2 



Rosenbluth 

re- analysis using Gep/Gmp from recoil polarization 



p{e,e') 



Gep/Gmp p{e, e'p) 



0.005 - 0.055 
0.13 - 1.75 
0.37 - 5.54 



Rosenbluth 
Rosenbluth 



|ll,[l6|,|37|,|l0T|| using Gmp from Ref. [|| 



G 



Mn 



d{e, e'n) 
d{e,e'N) 
d{e,e'N) 
d{e, e'N) 

d{e,e') 
d{e,e') 
^He{e, e') 



0.11 - 0.26 
0.12 - 0.61 
0.095, 0.126 
0.24 ~ 0.78 
0.07 - 0.89 
1.75 ~ 4.0 
2.5 - 10.0 
0.1, 0.2 



absolute, efficiency from d{'y,pn) 

ratio method, efficiency from p(7,7r+)n 

ratio method, efficiency from elastic p{n,p)n 

ratio method, efficiency from elastic p{n,p)n 

ratio method, efficiency from elastic p{n,p)n 

quasielastic 

quasielastic 

Fadeev analysis based upon Refs. [plJsl 



Geu. 



d{e, e'n) 
d{e, e'n) 
d(e, e'n) 
d{e, e'n) 
^He(e, e'n) 
^He(e, e'n) 

^20)^20 



0.26 

0.15, 0.34 

0.2 

0.5 

0.4 

0.67 

0.008 - 1.64 



m 



PWIA 

FSI analysis by Ref. |T03 | 
FSI from Arenhovel et al. 
FSI from Arenhovel et al. 



m 



Fadeev analysis by Ref. |55|| 
PWIA 

extracted from deuteron quadrupole form factor 



e(n, n) 







thermal neutron transmission in liquid 208 pj-, 
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TABLE II. Moments, radii, and per point for the four fitted densities. 



quantity 


model 


A 


xVn 


Mo 


M2 (fm2) 


R (fm) 


C (fm) 


Gev 

i-iy 


LGE 





0.67 


1.003 ±0.001 


0.776 ±0.019 


0.879 ±0.011 


0.879 ±0.011 






1 


0.69 


1.003 ±0.002 


0.712 ±0.020 


0.843 ±0.012 


0.881 ±0.012 






2 


0.71 


1.003 ± 0.002 


0.649 ± 0.021 


0.804 ± 0.013 


0.883 ±0.014 


Gev 


FBE 





0.69 


1.003 ± 0.001 


0.776 ± 0.020 


0.880 ±0.011 


0.880 ±0.012 






1 


0.70 


1.003 ± 0.002 


0.713 ± 0.024 


0.843 ± 0.014 


0.882 ±0.015 






2 


0.73 


1.003 ± 0.002 


0.651 ± 0.030 


0.806 ±0.019 


0.884 ± 0.020 


Gmp 


LGE 





0.71 


1.012 ±0.008 


0.729 ± 0.038 


0.849 ±0.023 


0.849 ± 0.025 






1 


0.71 


1.012 ±0.008 


0.659 ± 0.040 


0.807 ± 0.024 


0.847 ± 0.027 






2 


0.73 


1.013 ±0.008 


0.599 ± 0.037 


0.769 ± 0.024 


0.851 ± 0.026 


Gmp 


FBE 





0.71 


1.012 ±0.008 


0.729 ± 0.038 


0.849 ± 0.023 


0.849 ± 0.025 






1 


0.71 


1.012 ±0.008 


0.661 ± 0.040 


0.808 ± 0.025 


0.848 ± 0.027 






2 


0.75 


1.015 ±0.009 


0.615 ±0.073 


0.778 ± 0.046 


0.859 ±0.050 


Gmu 


LGE 





2.69 


1.016 ±0.025 


0.839 ± 0.099 


0.909 ± 0.055 


0.909 ± 0.058 






1 


2.69 


1.027 ±0.028 


0.823 ±0.116 


0.895 ± 0.064 


0.931 ± 0.067 






2 


2.71 


1.023 ± 0.027 


0.736 ± 0.107 


0.848 ± 0.063 


0.922 ± 0.065 


Gmu 


FBE 





2.69 


1.017 ±0.026 


0.843 ±0.101 


0.910 ±0.056 


0.910 ±0.059 






1 


2.70 


1.028 ± 0.029 


0.829 ±0.120 


0.898 ± 0.066 


0.934 ± 0.069 






2 


2.75 


1.026 ± 0.028 


0.752 ± 0.123 


0.856 ± 0.074 


0.930 ± 0.076 


Geu 


LGE 





0.52 




-0.115 ±0.003 










1 


0.55 




-0.115 ± 0.003 










2 


0.57 




-0.115 ±0.003 






Geu 


FBE 





0.55 




-0.115 ±0.003 










1 


0.58 




-0.115 ±0.003 










2 


0.56 




-0.114 ±0.003 
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